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Q Abstract 

We discuss the sensitivity reach of a neutrino factory measurement to non-standard neutrino 
interactions (NSI), which may exist as a low-energy manifestation of physics beyond the Standard 
' ^ ' Model. We use the muon appearance modes Ve — > ^n/i^e i^fi and consider two detectors, one 

at L = 3000 km and the other at L = 7000 km; The latter is nearly at the magic baseline which 
is known to have a great sensitivity to matter density determination. Assuming the effects of 
^ NSI at the production and the detection are negligible, we discuss the sensitivities to NSI and 

I— —I the simultaneous determination of ^13 and 6 by examining the effects in the neutrino propagation 

of various systems in which two NSI parameters Eap are switched on. The sensitivities to off- 
^ diagonal e's are found to be excellent up to small values of ^13. At sin^ 2^13 = 10~^, for example, 

{serl — a few X 10~^ at 3a CL for 2 degrees of freedom, whereas the ones for the diagonal e's are also 
^\ acceptable, leeelderrl) — 0.1(0.2) at the same CL. We demonstrate that the two-detector setting 

'""J is powerful enough to resolve the 013-NSI confusion problem, a notorious one which is thought to 

On be an obstacle in determining ^13 and 6. We believe that the results obtained in this paper open 

the door to the possibility of using neutrino factory as a discovery machine for NSI while keeping 
<^ its primary function of performing precision measurements of the lepton mixing parameters. 
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I. INTRODUCTION 



In the last 10 years various neutrino experiments made it clear that neutrinos have masses 
and lepton flavors mix [H, [2] . Though the measurement of the values of di^ and the CP phase 
6 as well as determining the neutrino mass hierarchy still elude us, there exist an array of 
experimental plans and ideas to make the goal. If successful they will bring us a more or less 
complete picture of how neutrinos organize their mass spectrum and how lepton flavors mix. 
When it happens, then, we will really be entering into the era of precision measurement of 
lepton mixing parameters. An intriguing question is whether this will serve just for precision 
measurements of more or less known quantities, or could help us to discover new physics 
beyond the Standard Model which is amended to incorporate neutrino masses and mixing. 

It has been suggested since long time ago that neutrinos may have non-standard in- 
teractions (NSI) in addition to those dictated by the Standard Model of electroweak in- 
teractions [21 in El E]- If there exists an energy scale for new physics at the TeV range, 
it is conceivable that it produces higher dimension operators which affect the way how 
neutrinos interact with matter, the point raised in [7]. It was argued that such NSI pos- 
sessed by neutrinos need not be subject to the stringent constraints that the charged lep- 
tons have to obey |8j. The existing constraints on NSI for neutrinos are worked out in 
[9l [To]. If we follow naive dimensional counting with four Fermi operators the order of 
magnitude of NSI is expected to be {mz iM^siY ~ (~ 10"'') for the energy scale 

of new physics Mjsisi of ~ 1 (10) TeV pjj. If we want to probe into such a tiny effect 
of NSI in the neutrino sector, extremely high sensitivities are required for experiments to 
detect deviations from the standard three flavor oscillations. A large number of papers 
in the literature is devoted to discuss NSI in various context including long-baseline ac- 
celerator experiments, the atmospheric and the solar neutrino observations, and supernova 
neutrinos [ni[l2l23l[ia[l5l[l6lEl[l8l[l9lEQl[2ll[22lEa 

It is the purpose of this paper to discuss the discovery potential of NSI in a neutrino 
factory. As is well known, a neutrino factory [27] is an ultimate apparatus for precision 
measurements of the lepton mixing parameters. Because of the capability of clean detection 
of muons including the charge identification, its sensitivity to is expected to go down to 
an extremely small value, somewhere in the range sin^ 2^13 = 10^^ — 10~^ [251 [SD]- Then, 
it is natural to think about neutrino factory as a discovery machine for NSI. 

One of the key issues in the discussion of the discovery reach of NSI is that we have to 
guarantee, at the same time, the ability of accurately measure the lepton mixing parameters, 
especially ^13 and 5, is left intact. However, it is a difficult goal to achieve in particular when 
we go down to very small values of and NSI contributions. In fact, it is known that the 
presence of NSI can confuse the measurement of by mimicking its effect [131 E] • Unless 
this problem is somehow solved, it is difficult to think about a neutrino factory as a sensitive 
hunting tool for NSI and at the same time as a viable apparatus for precision measurements 
of the lepton mixing parameters. 

We do overcome the confusion problem by our set-up. We consider a two detector setting 
with baselines at L = 3000 km and at L = 7000 km [T3j, which will sometimes be referred 
as the intermediate and the far detectors, respectively. Comparing the yields taken by the 
two detectors is the key to resolve the confusion. These features will be demonstrated in 
Sees. |v] and VI We must note, however, that it is only by keeping the solar Am^ terms that 
one can gain access to the CP phase 5 which is crucial to resolve the muddle between ^13 
and NSI. Therefore, we keep the solar Am^ terms throughout our analytical and numerical 
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treatment. 

If the two detector setting is a special requirement for hunting NSI it would not be easy to 
implement it into the project. Fortunately, it is not the case. It has been suggested |3Ti 
that combining measurements by detectors at two baselines, one at 3000-4000 km and the 
other at ~7000 km, is powerful enough to resolve the parameter degeneracy [311 E21 El] 
associated with the measurement of the lepton mixing parameters, a notorious obstacle to 
its precision determination. The power of the two detector setting has been confirmed also 
by more recent analysis [35]. 

In Sec.[TT| we present the general framework for probing NSI through neutrino propagation 
in matter by switching on two NSI contributions simultaneously and discuss some important 
points about NSI effects in the P(z/e —>■ z/^)/P(z/e —>■ v^) oscillation probabilities, the channels 



explored in this paper. In Sec. Ill, we elucidate the powerfulness of the strategy to look 
for NSI in a neutrino factory by combining the measurements of two identical detectors at 
different baselines. The assumptions we make about the neutrino factory parameters and 
the two detector setup as well as the analysis we perform are described in Sec. IV, In Sec. |V[ 
we discuss the maximal sensitivity to the various combinations of NSI parameters and their 
effects on the precision measurements of the oscillation parameters 5 and ^13. In Sec. VI we 



discuss the accuracy of the determination of the NSI parameters as well as their impact on 
the measurement of 5 and ^13 in the case they are sufficiently large to be measured by the 



neutrino factory experiment. The discovery reach to NSI, 5 and ^13 is presented in Sec. VII 



Sec. VIII is devoted to our final concluding remarks. In Appendix |A| we give details on the 
derivation of the appearance probability P(z/e z/^) in the presence of two NSI parameters. 



II. NON-STANDARD INTERACTIONS OF NEUTRINOS 



General features 



We consider NSI involving neutrinos of the type 



/-NSI 



(1) 



where Gp is the Fermi constant, and / stands for the index running over fermion species in 
the earth, f = e,u, d, in which we follow [9J for notation. P stands for a projection operator 
and is either = |(1 — 75) 01 Pr = ^(1 + 75)- We summarize here the bounds on 
which are obtained in P, [10] for the readers convenience: 



-0.9 < See < 0.75 ke^l ^. 3.8 X 10"^ 



-0.05 < e^^ < 0.08 \e 





< 0.25 


^^T 1 


< 0.25 




< 0.4 



(2) 



bounds from Davidson et al. (LEP) are at 90% (95%) CL. 

In this paper, we consider the effect of NSI in neutrino propagation in matter. It is known 
that NSI can affect production and detection processes of neutrinos so that a complete 
treatment must involve also the latter two effects. However, we should note that the muon 
storage ring which we assume as a source of neutrinos is special as it is one of the cleanest. 
Muon decay has been studied with great precision and room for NSI is smallest among 
the sources [36]. Because construction of the muon storage ring will require an intense 
muon source, we assume that the already tight constraints on NSI by muon decay would 
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become much more stringent at that time. Therefore, we beheve that to neglect NSI in the 
production of neutrinos in fact gives a fair approximation, unless we go down to extremely 
small values of 6^/3 • 

With regard to NSI at the detection this may be more debatable. However, we call 
the readers' attention to the fact that upon construction of the neutrino factory the near 
detector sitting in front of the storage ring will give a severe bound on NSI [9]. We also 
remark that even before that era several low energy neutrino experiments may be able to 
place equally severe constraints on NSI [371 EH EH]- The bounds from them are placed 
on the product of NSI at the source and the detection. But if the constraint on NSI by 
muon decays is strongly constrained, this will translate into a stringent bound on NSI at 
the detection. Furthermore, since we compare two identical detectors, at the intermediate 
(L = 3000 km) and the far (L = 7000 km) location, the NSI effect at detectors tend to 
cancel. Therefore, our approximation of ignoring NSI at the detection points may give a 
reasonable first approximation. 

To discuss the effects of NSI during neutrino propagation in matter we will use the effec- 
tive coefficients Sa/s as it is traditional in this field. They are defined as Eafs = ^/p ^^a^- 
where n/ is the number density of the fermion species / in matter. Normalizing by 
the electron number density He leads to a simple structure of the effective Hamilto- 
nian which governs the neutrino propagation in matter. Approximately, the relation 
Ea/B — + 3 E^p + 3 holds because of a factor of ~3 larger number of u and 

d quarks than electrons in iso-singlet matter. 

The time evolution of neutrinos in flavor basis with non-standard neutrino matter inter- 
actions can be generically written as 
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(3) 



where U is the Maki-Nakagawa-Sakata (MNS) [Tj matrix, and a = l^'^GpUf.E ^ where 
E is the neutrino energy and denotes the electron number density along the neutrino 
trajectory in the earth. Amf^- = — with neutrino mass mj = 1 — 3). Eq. (j3]) defines 
the framework for discussing neutrino propagation in matter with NSI. 



B. Physics of neutrino propagation in matter with NSI; Two e system 

In this paper, we will be dealing with the cases where the following pairs of e parameters 
are present at the same timei e^-y e^^^ ^er ^tt-i ^tt ^ee? ^ee? ^tt-j and e^j- e^^. 
The systems with many e's are complicated enough and a step by step approach is needed 
to grasp the whole picture. In Appendix |X| we will derive the exact expressions and the 
approximate formulas for the appearance oscillation probabilities P(z/e — >■ z^^) in the two 
systems Eer — £ee and Eg^ — E^e- To obtain tractable formulas we use perturbation expansion 
in terms of small parameters, sis ^ ^21 — ^er (or Ee^}j which we assume to have comparable 
sizes ~ 10^^, where is defined in Appendix [Aj For simplicity, we collectively denote their 
order of magnitude as e under the hope that no confusion arises with the NSI elements E^ji- 
See [HI Uni EHl HO] for different treatment of the perturbation series. We will observe 
that in leading order the appearance oscillation probabilities Pive ^fi) are of order e^. Of 
course, it reduces to the Cervera et al. formula [28] if 5^/3 are switched off. 
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We mention some notable features of the leading order formulas of P(z/e v^) in (A25) 



and (A27), and make comments for better understanding of our results. To discuss simulta- 
neously the Ser — See and Se^ — See systcms, we use the notation for (a, /3) is either (e/x) or 
(er). Since an off-diagonal element of e can have CP violating phase we have a new source 
for CP violation in addition to 5 in the MNS matrix, as emphasized in [TT] in a somewhat 
different framework. In this context we notice that Sap comes in via three different ways: 

(a) SafB comes into the formula with its phase, 0q,/3, in the form 5 + (j)ai3 together with the 
leptonic Kobayashi-Maskawa (KM) phase S in the MNS matrix. fTEl . 

(b) ea/3 comes in by itself, i.e., in the form of Re^Eap) or Im^Eap)- 

(c) Eais comes in by the absolute magnitude squared, |ea/3p without phase. 

Terms of the type (a) and (c) survive even when the solar Am2i is turned off, whereas terms 
of the type (b) arise only when they are accompanied by the solar Amgi. 

We suggest an intuitive understanding of the fact that two phases come together when 
Am2i is switched off. The system without Am2i contains effectively only two generations 
of neutrinos and on physics ground the CP violating phase should be unique. Therefore, 
the phase of Eap and the KM phase 6 must come together. From this reasoning, we suspect 
that this feature holds even in the exact formula. Of course, the two phases start to play 
their individual roles when the solar /S.m\^ is turned on, as in the form of (b) above and the 



second line of the Cervera et al. formula (A26), as seen in (A25) and (A27). 



Due to the special way through which the phase of Eafi enters into the oscillation proba- 
bility in the absence of the solar Am2i, we expect that confusions occur in the experimental 
settings for which ignoring the solar Amgi gives a fair approximation. We expect two types 
of confusion to occur: 

• Two phase confusion: Since they appear in the special combination (^^p + ^i it is 
obvious that two phase confusion occurs; the data allows determination of only the 
sum. The confusion cannot be resolved even if the anti-neutrino channel is combined. 
It can be resolved for relatively large \Eai3\ because the terms with /^m^^ start to play 
more important role. These features are in fact observed for Ee^j, in [25j . We note that 
the symmetry under simultaneous shift of and 5 by the same amount but with 
opposite signs is a symmetry of the whole system at the magic baseline. We will see 
that this confusion shows up in our analysis. (See Sees. |V]and VI ) 



Phase-magnitude confusion: We expect that another type of confusion exists. The 



term which contains the phases in (A25) and (A27) take the form 



A\Eal3\ COs(0„/3 + (5) + B\Eap\ sin(0a/3 + 5) = \EapW + 5^ COs(0a^ + ^ - > (4) 

where tan^ = B /A. Therefore, varying the magnitude of the NSI element can be com- 
pensated by adjusting the phase of Eap- Unlike the two phase confusion, this confusion 
can in principle be resolved if (1) |£a/3p terms are important enough to resolve the con- 
fusion, and/or (2) the anti-neutrino channel is combined, because the coefficients A 
and B are different between neutrino and anti-neutrino oscillation probabilities. 
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III. PROBING NSI BY DETECTORS AT TWO BASELINES; THEIR CHARAC- 
TERISTICS AND SYNERGY 



In this section, we indicate by qualitative level arguments that the v^jv^ appearance 
measurement in neutrino factory at the magic baseline, 

and the synergy obtained when combined with measurement at L ~ 3000 km provide a 
powerful method for probing non-standard neutrino-matter interactions. Powerfulness of 
the measurement at the magic baseline may be natural to expect for e^e because it is effec- 
tively equivalent to measuring the electron number density in earth matter for which the 
accuracy of determination is known to be excellent [HI |12]. We show in this paper that this 
statement is even more true for flavor off-diagonal NSI. We rely on v^. v^jve appear- 
ance channels, the so called golden channels [.28j, because of the matured muon detection 
technology. 



A. Detector at the magic baseline as a sensitive probe to NSI 

To illuminate that measurement at the magic baseline is a powerful indicator for NSI, in 
particular for off-diagonal Eap, we present in Fig. [Tjthe ellipses formed when 6 is varied in 
the bi-probability space spanned by P(z/e —>■ v^) and P{i'e ^ii) [SS| for various non-zero 
NSI parameters. The calculation was performed numerically assuming the constant matter 
density p = 4.5 g/cm^. In each panel only the indicated particular Sais is turned on. Also 
shown in Fig. [l] as an orange strip is the region covered by the ellipses when ^13 is varied 
for the cases without NSI. It approximately forms a narrow straight strip (or "pencil") 
because of vanishingly small effect of S at the magic baseline |13]. The remarkable feature 
of Fig. [1] is that the effect of NSI in the neutrino and the anti-neutrino probabilities is large 
for electron-type off-diagonal terms, Eer and Ee^, even though their size is extremely small, 
p „ = e^-r = 10"^. Notice that the sizes of the other e's which give rise to effect of similar 
magnitude as those of Ee^ and Eer are larger by a factor of 20 (e^,- and e^^) and of > 50 
(Err)- This is the key point of our setting which allows the extremely high sensitivity to e^t 
and 

There are some curious features in Fig. [l| The behavior of the ellipses with Eer and Se/i 
are distinct from the other cases having ellipses shrunk forming almost lines. At the same 
time it is also notable that they look almost identical to each other. Let us understand these 
characteristic features of Fig. [T] 

In Appendix [A] we derive the leading order formula for P(z/e v^) with Eer 

or S(>^ cis well 

as Eee corresponding to the ones derived by Cervera et al. [28j for the standard case without 
E. At the magic baseline, f§ = tt, the formulas for the appearance oscillation probability 
greatly simplifies. With NSI represented by Eer it is given by 

Pip, -^ p„■E,^ = A-^ — H^^4rTTTslsf,sm^ ' 



^ 4ac23si3 



a — AmU^ 



2/\ml-ySi^\Eer\ C0S(5 + (per) + C23akeTp siu^ f ^'^ \ . (6) 
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IJi-Pi-iibjljilily pkL Joi- Lhc Magic IJasclmc l^T2DQkrn, E=30 GeV fnr sin 2e^^ = 0.CCO5, 0.00 1 :inii O.OOl.'i 



0.20 



^0.10 



I I I I I r 1 f 1 I I I 

le 1=0.02 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
■IE 1=0,001 ■ 




0.1 0.2 OJ 0.4 0.5 



0.1 0.2 0.3 0.4 0.5 



0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 1: Bi-probability plots in P(fe — > f^) — P{i^e — > t'/x) space at the magic baseline, L = 7200 
km, for E = 30 GeV, computed numerically using the constant matter density p = 4.5 g/cm^ with 
the electron number density per nucleon equals to 0.5. The both axes is labeled in units of 10~^. 
In each panel only the indicated particular e^^ is turned on. The upper (lower) panels, from left to 



right, correspond to the case of and e 



respectively. The red 



and the blue ellipses are for positive and negative signs of e, respectively, for the cases with (from 
left to right) sin^ 2^13 = 0.0005, 0.001, and 0.0015, as indicated in the heading. The values of the 
non-vanishing e are written in each panel. The orange colored region indicates the region spanned 
by ellipses without NSI when ^13 is varied. The green dots are unresolved ellipses corresponding to 
the same values of sin^ 2^13 but without NSI. The values of the standard lepton mixing parameters 
are given in the caption of Fig. [sj Only for the case of e > and sin^ 2^13 = 0.001 we show the 
position corresponding to the four different values 6 = 0,7r/2,7r and 37r/2 by the open circle, 
square, diamond and asterisk, respectively. 



The corresponding formula for anti-neutrinos can be obtained by making the replacement 
a —a, 5 — * —S, and (per —(per- The formula with is very similar to ([6]). It is 
obtained by replacing C23£eT by S23£efi in the second line of Eq. Then, the value of 
P(z/e z/^) with Est and Eefj. are numerically equal with the values of the parameters used 
in Fig. [l| the maximal value of ^23 and Iser] = ke/x|- It explains the identical behavior of 
the probabilities with and Eer seen in the second and the third upper panels in Fig. [T| 
The ellipse in these panels shrink approximately to a line because there is only cos (5 + ^er) 
dependence in the probability where 0er denotes the phase of Eer- The form cos((5 + (per) 
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BL-Pn)babLUby plot for the boBclinc L=3000km, E=30 GcV for sin'^S^^ = 0.0005, O.OOl and 0.CD15 




0.5 1 1,5 2 0,5 1 1.5 2 05 1 1,5 2 2.5 

P(v-v^) [jtlO *] 

FIG. 2: The same as in Fig. [T]but for the baseline L = 3000 km with the matter density p = 3.6 
g/cm'^. The same values of e are used in each panel. In the left and right lower panels the ellipses 
with positive and negative sign of e overlap almost completely and each individual curve is not 
visible. The green ellipses which correspond to the same three values of sin^ 2^13 but without NSI 
are clearly visible. 



means that the shape and the location of shrunk ellipses in Fig. [T] are unchanged even when 
(per is varied. The label of 5 on the ellipse (if it is placed), of course, changes. 

The length of the shrunk ellipse is equal to twice the coefficient of cos((5 + ^er) of the 
last term in Eq. ([6]). With the values of the parameters used (E=30 GeV, L = 7200 km, 
p = 4.5 g/cm^, Egr = 0.001, sin^ 26*13 = 0.001) it can be estimated as 1.8 x 10^^, which is in 
perfect agreement with the length of the shrunk middle ellipse projected onto P(i^e '^fi) 
axis in Fig. [TJ Therefore, we have understood the qualitative features of the flat ellipse at 
the magic baseline as well as its size. 

We can also understand the reason why the effect of Err and are suppressed by deriving 
a similar formula which contains one of them in an analogous way as in Appendix |X] The 
leading order terms which involve one of these e's are of order and hence smaller. 
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B. Detector at L = 3000 km and the synergy expected when combined with the 
one at L ~ 7000 km 



It is instructive to compare the similar plot for the intermediate detector at L = 3000 
km which is shown in Fig. |2l to Fig. [T] at L = 7200 km. The same values of e's as in each 
corresponding panel in Fig.]l]are used. The three representative values of sin^ 2^13 used are 
also the same as in Fig. [T} We immediately notice several clear differences. First of all, the 
effect of 6 is large at L = 3000 km for both cases with and without NSI, as indicated by 
the blue and the red lines (with NSI) and by the green ellipses as well as by wide span of 
the ellipses indicated by the orange region (without NSI). The ellipses with NSI, except for 
the case with Sg^, are buried into this region. This is clearly the cause of the problem of 
confusion between ^13 and e's one encounters when one tries to measure ^13 and S allowing 
for NSI; The system with e's can mimic the one with different values of ^13 but without NSI. 
In Fig. |2] the difference between the confused parameters and the genuine ones are small, 
apart from the because we take small values of e's. It is also notable that the effect 

of the sign of the e's is not quite visible, which is nothing but a consequence of the two-phase 



degeneracy as explained in Sec. II B (See Sec. VI B for more about the degeneracy.) 



When we try to determine the values of e's it is also a bad news because there can be 
a severe confusion between NSI and the standard oscillation effect with 613 and 6. We will 
see in Sees. |V] and IVII that this confusion for e determination is much severer than that in 
^13 — S determination at L = 3000 km, which will be manifested as a complicated island 
structure in the allowed region. On the other hand, when we try to measure ^13 and S at 
L = 3000 km confusion due to the presence of e's is not so significant provided that the e is 
small as we see in Fig. [2j In particular, as we will see, the sensitivity to 6 is good even after 
marginalization of e's. (See Fig. [7]-[To]in Sec. |Vj) 

On the other hand, at L ~ 7000 km the sensitivity to e's is great though essentially there 
is no sensitivity to 6. These consideration naturally suggest the possibility of combining the 
intermediate and the far detectors to determine simultaneously NSI parameters at the same 
time measuring ^13 and 6. 

The markedly different behavior of the bi-probability plots between the e^^ and sys- 
tems at L ~ 3000 km can be traced back to the difference between the third terms in the 
analytic formulas (A25) and (A27). We should note that the almost identical behavior be- 
tween the systems with Ee/i and Ser sA L = 7000 km, and the marked difference between 
them at L ~ 3000 km makes it interesting by itself to compare their sensitivities at these 
two baselines and when they are combined. We will see later that sensitivity to Eefj. is es- 
sentially determined by the effect at L = 3000 km whereas that to is determined by the 
combination of L = 3000 km and 7000 km. 



IV. ANALYSIS METHOD 
A. Assumptions 

In our analysis in this paper, we make the following assumptions for the parameters of 
neutrino factory. We assume an intense muon storage ring which can deliver 10^^ useful 
decaying muons per year. The muon energy is taken to be 50 GeV. We assume 4 years 
running in neutrino and 4 years running in anti-neutrino modes, respectively. We assume 
the two magnetized iron detectors, one at baseline L = 3000 km and the other at L = 7000 
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km which is close to but not exactly the magic baseline Lmagic-^ Each detector is assumed to 
have fiducial mass of 50 kton. We consider the golden channels, Vf> — > i/^ and Vf. — >■ z/^ in this 
paper. For simplicity, we use the constant matter density approximation throughout this 
work, and take the Earth matter density along the neutrino trajectory as p = 3.6 g/cm^ and 
p = 4.5 g/cm^ for baselines L = 3000 km and at L = 7000 km, respectively. The electron 
fraction Ye is assumed to be 0.5. We believe that using more realistic Earth matter density 
profile will not change much our results. 

In most part of our analysis we make the following simplifications: (1) We ignore all 
the background due to the misidentification of muon charges and other causes. (2) We 
neglect the systematic uncertainties. Since the muon detection at high energies is supposed 
to be extremely clean in magnetized iron detectors the simplification (1) may not affect the 
results in a significant way apart from the case of extremely small di^. With regard to the 
systematic errors we feel that no solid numbers are known yet in spite of the fact that great 
amount of efforts are made toward reliable estimation of them [H]. Nevertheless, we will 
try to estimate to what extent the sensitivities we obtain in our analysis are affected by 



introduction of background and the systematic errors. See Sec. VE 



We always take the normal mass hierarchy as an input. We will consider in this paper 
various cases in which only two different flavor elements of NSI (^a/j) are turned on at the 
same time. We also assume that 5^/3 is real, leaving the interesting topics of complex phase 
effects to elsewhere. 



B. Analysis procedure 

We deflne the function as follows, 

where N°j^^ is the number of observed (simulated) events computed by using the given input 
parameters and N^^k is the theoretically expected number of events to be varied in the flt by 
freely varying the mixing and NSI parameters. Since we ignore the systematic uncertainties, 
the denominator in ([T]) represent the statistical uncertainties. The summation with respect 
to indices i,j and k imply energy (3 bins), baseline (3000 km or 7000 km), and the type 
of neutrinos (neutrino or anti-neutrino), respectively. The intervals of 3 energy bins we 
consider are 4-8 GeV, 8-20 GeV and 20-50 GeV for neutrinos and 4-15 GeV, 15-25 GeV and 
25-50 GeV for anti-neutrinos. 



By assuming the far detector at baseline somewhat off Lmagic, we want to demonstrate that it is not quite 
necessary to place it exactly at Lmagic- First of all, we feel it unrealistic that one can place the far detector 
at Lmagic with a mathematical precision. Apart from the problem of site availability, the exact magic 
baseline cannot be determined prior measurement unless the values of the relevant mixing parameters, 
the earth matter density along the neutrino trajectory and its relationship with the effective density for 
neutrino oscillation |42j are precisely known. Only a posteriori the matter density can be measured in 
situ in the experiment |41j. 
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The theoretically expected number of events are computed as 

N''^^°{e,„6,e) =n^TM^^-^ / g{E) cr,^^,^){E) P,^_.,^^,^^^,^){E;d,„6,6)dE, (8) 

where is the number of useful muon decays per year, T is the exposure period (in years), 
M is the detector mass (in ktons), Na is the Avogadro's number, is the muon mass, Efj_ is 
the energy of the stored muons, L is the baseline, a^^(p^){E) is the charged current interaction 
cross section for i/^ and z/^, and Pi,^^^^(^p^^p^){E; 613, 6, e) is the oscillation probability. In 
this work, we considered the case where E^^ = 50 GeV, M = 50 kton, T = 4 yr for both 
neutrinos and anti- neutrinos, and = 10^^ per year. The function g{E) which is given as 

is the unoscillated z/g or z/g energy spectrum normalized to 1. 

We assume, for simplicity, the detection efficiency is 100%.^ We neglect the finite energy 
resolution in the detectors. Since the number of energy bins are small (=3) inclusion of the 
energy resolution will not alter the results in a significant way. 

The observed number of events are computed exactly in the same way but using the given 
input parameters of ^13, S and NSI parameters (e), so that Xmm = at the best fit point by 
construction. 

Using the function defined in Eq. ([7|, we define the allowed (sensitivity) regions by the 
commonly used condition, Ax^ = ~Xmm = 2.3, 6.18 and 11.83 for 1, 2 and 3 a confidence 
level (CL) for 2 degrees of freedom (DOF), unless otherwise stated. 



V. SENSITIVITY TO NON-STANDARD INTERACTIONS AND MEASURE- 
MENT OF AND 6 WITH NSI 



In this section, we discuss the sensitivity to NSI by the two detector setting. We also 
discuss the accuracy of the determination of the mixing parameters ^13 and S in the presence 
of NSI. For these purposes we take the input values of two e's to be zero (or equivalently 
vanishingly small) but freely vary them in fitting the data. The case with non-zero input 
values of e's will be dealt with in Sec. VI To demonstrate the synergy between the interme- 
diate and the far detectors we present the sensitivity to NSI for each detector as well as the 
combined one throughout this and the next sections. As a typical value of di^ we consider 
the case of sin^ 26^13 = 10~^, though we also discuss the case with sin^ 26'i3 = 10""^ to show 
how the sensitivities depend on ^13. We have examined the four values of 6, n/A, 7i/2, vr, and 
37r/2 as representative cases. However, we present only part of the figures we have drawn 
not to make this paper too long. 



^ If the efficiency is /xIOO% we are effectively assuming the fiducial mass of the detector of 50// kton. In 
the current estimate / is expected to be about 0.8 apart from an extremely low energy region ~ 5 GeV 
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A. Constraining NSI; Case of zero input 



In Figs. [3]j6j we present the constraint on NSI that can be imposed by the neutrino factory 
measurement defined in Sec. IV A| with various selected combination of two e parameters.^ 
Fig. |3] and |4] are for 5 = 7r/4 whereas Fig. |5]and|6]are for 5 = 37r/2. The left, the middle, and 
the right panels of Figs, [s] and [s] (Figs. |4] and |6]) are for the combination, showing horizontal 
- vertical axes, E^e ^eri ^tt ^er; ^ee ^tt (^ee ^e/x; ^tt ^e/iy ^er ^efi): respectively. The 
top, the middle, and the bottom panels are, respectively, for the baselines L = 3000 km, 
L = 7000 km, and the two baselines combined. The blue, the red, and the green curves are 
the allowed contours at Icr, 2a and 3cr CL for 2 DOF, respectively. 

We observe that from the results presented in Fig. [3] through Fig. |6j the detector at 
L = 3000 km alone does not have good resolution power for the possible existence of e's, 
except for eg/^- The parameter ^e^t is special since its impact on the oscillation probability is 
so large (as seen in Fig. [2]) that the measurement at L = 3000 km alone can give a strong 
constraint on which seems consistent with the result obtained in Ref. [25]. Apart from 
the cases which involve Eefj., the effect of the simultaneous presence of two e's is manifest in 
the appearance of many correlated regions/ islands, though the precise shapes of the regions 
depend on which combination of e's is turned on and on which input value of 6 is used. 

A similar statement applies to the case of the detector at L=7000 km. But, the correlation 
between e's is quasi one-dimensional for most of the combinations though we see branch 
structure in Eee — Ser case and sizable width in the Srr Overall, the constraints 

on the diagonal elements. See and Err, are much looser compared to those on e^r and 6^^. 
The latter feature comes from high sensitivities to the off-diagonal e's expected at the magic 
baseline as discussed in Sec. Ill A while less sensitivity to the former is expected by Figs. [T] 
and [2j The feature of correlation is least obvious in the combination See — Stt though 
indicating a weak oblique linear dependence. For the branch-like structure we will make 
comments to clarify its nature in Sec. VI because the structure is even more prominent with 
non-zero e input. 

The effect of combining the intermediate and the far detectors is remarkable. The allowed 
regions scattered in wide ranges in the top (3000 km) and the middle (7000 km) panels 
combine into a much smaller region in the bottom panel in Figs. |3]j6j We should remark 
that although the above mentioned over-all features remain unchanged for different values 
of 5, the resultant sensitivity to e's depends rather strongly on the value of CP phase 5 as 
one can see by comparing between Fig. [3] and Fig. |5| The problem of 5 dependence of the 
sensitivity to e's will be fully addressed in Sec. |VII[ 

One may ask the question why the scattered regions in the top panel and the extended 
region in the middle panel in each column in Fig. |3] and Fig. |5j and Fig. |4] and Fig. |6j 
can be combined to yield such a small region. The answer is that it is due to the CP 
phase 5. Namely, most of the region of overlap between the top and the middle panels 
have mismatch in value of 5, and hence they do not survive when the two constraints are 
combined. Therefore, keeping the solar Am? and the KM phase degree of freedom is the key 
to the high sensitivity to e's we observe in Figs. |3]j6| The synergy effect that merit us by 
combining the intermediate and the far detectors are even more significant compared with 



^ It appears to us that ragged behavior of the contours seen in some of the plots is of physical origin due 
to the complicated structure of four-dimensional function. Some of the structures, however, could be 
smoothen to a certain extent by introducing a finer grid and the finite energy resolution. 
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the one in identical two detector method for measuring CP violation and determining the 
mass hierarchy [l5l |l6] . 

Our last comment in this subsection is that a characteristic feature that manifest itself in 
Figs. [3]j6] gives us a warning. Namely, one could significantly overestimate the sensitivity to 
detection of non- vanishing NSI by working only with a particular single element See or ^er, 
for example. By working with two e's at the same time one is able to recognize the whole 
structure as given in these figures.^ 



B. Sensitivity to 6*13 and 6 in the presence of NSI; Case of zero input 



We now examine the sensitivities to ^13 and 5 in the presence of NSI. In Fig. [7] through 
Fig. 10 we present the sensitivities to ^13 and 5 by taking the same input value of ^13, 
sin^26'i3 = 10~^. Figs, [t] and [s] are for 5 = 7r/4, whereas Figs, [ojand 10 are 6 = 37r/2. The 
organization of the figures is the same as in Figs. [3]j6l 



Let us look at the top panel for the detector at L = 3000 km. We clearly observe the 
phenomenon of "confusion" in the presence of NSI, in particular for the case of 6 = 37r/2. 
Namely, NSI can mimic the effect of non-zero ^13 so that the allowed region extends to a 
very small value of ^13. However, this feature is strongly perturbed by the measurement at 
L = 7000 km. Because it is highly sensitive to the effects of NSI, it helps to resolve the 
confusion between NSI and ^^13. After combining informations of the intermediate and the 
far detectors a tiny region in sin^ 26'i3 — 6 space results, as one can observe in the bottom 
panels in Fig. [7] through Fig. 10 Thus, the problem of NSI - 6is confusion can be solved by 
the two detector setting.^ 



C. How do the sensitivities depend on ^13? 



We discuss briefly how the sensitivities change when we take a smaller value of 6*13, 
for example, sin^2^i3 = 0.0001. We present our results only for 6 = 37r/2 for reasons 



of limitation of space but make some comments on the other cases. In Figs. 11 and 12 



we show the sensitivities to NSI parameters which may be compared to the corresponding 
figures. Figs. [5]and[6} for the case sin^ 26*13 = 0.001. To our surprise, the over-all features of 
the allowed regions are similar between the two cases. The only notable changes are that: 
(1) in the See — £tt system the sensitivities to these e's becomes worse by a factor of 3-4, 
but (2) in all the three systems which involve Eefj,, See — £eii, £tt — £efj., and Eer — £e^l, the 
sensitivities become better for sin^ 26*13 = 0.0001. The similar worsening or improvement 
of the sensitivities to NSI are observed in other values of 6. They become worse in some 
cases, in particular, in the See ~ Ett system with 6 = tt/A and 7r/2 at L = 3000 km. The 
sensitivities to See and Err improve by almost a factor of 2 at 5 = vr with sin^ 26*13 = 0.0001 



It also raises the question whether the two e systems are sufhciently generic to reveal the full structure 
of the systems with larger number of e's which have a multi-dimensional manifold. We are not able to 
answer the question in this paper. 
^ The "confusion theorem" stated in [T3] refers to the situation in which we have both NSI at source and 
in propagation and they are related with each other in a specific way. Therefore, our discussion in this 
paper does not affect the validity of this theorem. 
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compared to the case of 0.001. The sensitivity to Ser gets improved by about a factor of 2 
in the ^e/^ ~ ^er system at 5 = ir/A. Thus, most of the improvement occur in the systems 
which involve £e/i- The reason is that the sensitivity to ^e^i is particularly good at L = 3000 
km and it makes the synergy effect even stronger. 



We notice that the branch structure we saw in Sec. V A is also rounded off, and the island 
structure becomes less prominent at sin^ 2^13 = 0.0001. The smoothing of the contour and 
the improvement of the sensitivity may be due to the fact that at such very small 61^ the 
effect of NSI becomes important compared to the standard oscillation effect. 

On the other hand, the sensitivity to ^13 and 6 definitely becomes worse when we go 



down to sin^ 2^13 = 0.0001 as shown in Figs. |13|and|14[ First of all, the problem of 6'13-NSI 



confusion becomes severer as one can see by comparing Figs. 13 and 14 to Figs. |9]and 10 A 
new allowed region emerges in the e^e — Ser and Srr — ^er systems at L = 3000 km. Due to 
lack of statistics at the small 6'i3 sensitivity to ^13 is significantly reduced at L = 7000 km, 
leading to the visible loss of the sensitivities. 

In summary, neutrino factory with two detector setup at L = 3000 and 7000 km is more 
resistant to small values of ^13 as a discovery machine for NSI, rather than as a machine for 
precision measurement of 6*13 and b. 

D. Comparing the sensitivities to ^13 and b for cases with and without NSI 

It may be worthwhile to compare the sensitivities to 6*13 and 5 for cases without NSI to 
the one with NSI using the same machinery used in our foregoing analysis. In Fig. 15 we 



present the results without NSI. The results are to be compared with those in Figs. [7 10 



The resultant sensitivities to ^^13 and b are extremely good compared to those obtained for 



the cases with NSI. We will quantify this comparison in a more extensive way in Sec. VII 
where we will cover the whole space of interest in ^13 and 5. 



E. Effect of background and systematic errors 

To have a rough idea of the possible effects of the background and the systematic uncer- 
tainties to the sensitivities we have repeated the computation with these effects at 5 = 7r/4 
for various values of 6*13 in ege-^er, ^rr-^eT-, and See-^TT systems, which correspond to Fig. |3| 
Fig. [7| and 5 = 7r/4 counterpart of Fig. 11 (not shown). We estimate the number of back- 



ground events by using the signal to noise ratio calculated in Ref. [47j, and rescaling the 
background to the number of useful muon decays used in this paper. We assign 20% un- 
certainty to the number of background events. We also take 2.5% as the total systematic 
uncertainty on the measurement. 

For sin^ 26'i3 = 0.001, we find that the introduction of systematic uncertainties and back- 
ground lead to decrease in the sensitivity to ^er by ~ 30%. The worsening of the sensitivity 
to Bee and depends upon the sub-system being within the ranges 30%-60% and 40%- 
50%, respectively. The effect of systematic errors and background is less prominent at 
sin^ 2^13 = 0.005 but becomes larger at sin^ 2^13 = 0.0001, making the sensitivity a factor of 
2 worse. 

On the other hand, the sensitivity to b and sin^ 26'i3 would be worsen by 30%-40% and 
50%, respectivelly, at sin^ 26*13 = 0.001, more or less independently of the sub-system treated 
and on the e's discussed. An interesting feature of the sensitivity to 5 and 6*13 is that the 
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effect of systematic errors diminishes as ^13 becomes small, and is at ~ 20% (30%) level for 
S and sin^ 2^13, respectivelly, at sin^ 26'i3 = 0.0001. Therefore, the impact of including the 
systematic to the sensitivity of 6 and ^13 is limited in size, which is good news. All numbers 
were evaluated at 3 cr CL. 

Though we continue to ignore the systematic uncertainties and the background in the 
rest of the analysis in this paper, the reader should remind that the estimated sensitivities 
to NSI, 6, and 6*13 have uncertainties at the level quoted above. To completely remove the 
uncertainties requires precise knowledges of the performance of the detector. 



VI. ACCURACIES OF DETERMINATION OF NSI, 613, AND 5 

We now discuss the question of how well the magnitude of non-vanishing NSI can be de- 
termined by the two detector setting, and at the same time to what extent the measurement 
of 613 and 6 can be affected by non-zero input values of NSI. In this section, we examine the 
particular systems: Eee—^er, ^TT — ^er, and Eee—^TT- We do not discuss the case which involve 
Ee^- The reason is that we want to take the NSI input values well below the current bounds 
presented Eq. ^ . The input value of e^^ which is comparable to the present constraint will 



not affect the sensitivity. In Sec. VII, we will present more global informations. We work 
with the particular input values of 613, 6, and the e parameters; sin^26'i3 = 10"'^, 6 = 3tt/2, 
Est = 0.01, Eee = 0.1, and Err = 0.2. But, the features are not so different for other input 
values. 



A. Sensitivity to NSI, Oi^, and 5; Case of non-zero input of NSI 

Roughly speaking, the sensitivities to NSI, 6^13, and 6 are not affected so much by the 
non-zero input values of e's. This is because, even in the case of vanishing input, we freely 
vary them in fitting the data, and the presence of these extra degrees of freedom of varying 
over the NSI parameters is of key importance to determine (or affect) the sensitivities. Since 
the results are similar to the previous case, where NSI input were set to zero, we only present 
the figures for the case S = 37r/2. Generally speaking, the sensitivity is worse than for the 
case 6 = 7r/4. More or less one can guess what would be the general characteristics of this 
case from the figures with zero input. Figs. |3} [4], [7|and[8j 



We first discuss the accuracy of the determination of NSI parameters. Presented in Fig. 16 
is the sensitivities to NSI parameters with the non-zero inputs. We notice that the features 
of the allowed regions are essentially the same as in the zero input case shown in Fig. |5| 
(About the two split islands in the middle panels see below.) We note that the accuracy 
of the determination of the e parameters varies case by case. Most significantly, the size of 
allowed region of Eer, £ee, and e^-t become worse by more than a factor 2 in Eee — ^er and 
e^-^ — Eer systems. However, the accuracy of the determination of Eer and Err at 6 = tt/A 
(not shown) improves by a factor of ~ 2 with the non-zero input values of NSI. 

Next, we discuss the accuracy of the measurement of ^^13 and 6 in the presence of non- 



vanishing input values of NSI. Fig. 17 serves for this purpose. Again the features of Fig. 17 
are similar to those of Fig. |9j We, however, notice some new characteristics. The ^13-NSI 
confusion at L = 3000 km is severer for 6 = 37r/2, but it is milder for 5 = 7r/4 (not shown) 
in the present case compared to zero e input. Therefore, the degree of confusion depends 
very much on 6. Also the resultant accuracy of the determination of 613 and 6, after the 
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two detectors are combined, depends upon which e is tuned on and also on 5 though only 
mildly. For 5 = 371/2, the uncertainty on sin^ 2^13 is smaller (larger) for the case of non-zero 
input value of the NSI parameters than for the zero-input case in the Err — ^er {^ee — ^er) 
system. It is interesting (and encouraging) to observe that non-zero input values of NSI 
essentially do not disturb the sensitivities to the determination of the NSI parameters but 
also the sensitivities to ^^13 and 5. 



B. Parameter degeneracies with and without NSI 



In this subsection we want to make some remarks on the degeneracies associated with 



the measurement of NSI parameters as well as 613 and 6. In Fig. [16} we observe two clearly 
separated islands in the regions allowed by the far detector measurement (middle panels) for 
the Err — £eT and the Eee — Ser systems. A similar structure exists also in the combinations 
involving Ee^- There exits a degeneracy of solutions of e's. It is easy to understand the cause 
of the degeneracy; It is due to the invariance under sign change of Est that can be absorbed 
by rotation of 5 by tt. The fact that the symmetry exists only at the magic baseline can be 
easily seen by comparing Eqs. (A25) a nd ([6 ). This is also true in the corresponding formula 



with Eeit that one can derive from Eq. (A27). Since 6 is marginalized the two clones appear 



in the allowed regions in the Err — ^er and 

When the measurement at the two detectors, intermediate and far, are combined we 
obtain a unique allowed region; The degeneracy is resolved. It must be the case because 



there is no such symmetry in the whole system, see Eq. (A25) or (A27). But, it is nice to 



see that it actually occurs in the concrete two detector setting adopted in this paper. 

We want to make brief remarks on the conventional parameter degeneracy [311 
which would affect the sensitivity to the lepton mixing parameter. It might also affect 
the sensitivity to NSI through its effect to ^13 and 6. First of all, there is no 623 octant 
degeneracy because 623 is held fixed to 7r/4 in our analysis. With regard to the four-fold 
^13 — 6 degeneracy duplicated by the unknown sign of Amg^^ we have no hint for its existence 
in our analysis, apart from some limited cases at sin^ 2^13 = 10~^ which might be affected by 
it. The reason for this observation is that the clone solutions have a significantly different 
value of S from the true solution (very roughly speaking 5cione ~ vr — 5true), in both the 
intrinsic and the sign-Amgj^ degeneracies. Nonetheless, from the figures we presented (and 
also in all those not shown) these is no such clone region. The reasons are: (1) The intrinsic 
degeneracy is resolved by the spectrum information^, and (2) With long enough baselines 
the mass hierarchy is determined. Therefore, we suspect that the conventional parameter 
degeneracy plays minor role, if any, in our sensitivity analysis, except for the one associated 
with 6*23 which we do not consider in this work. 



VII. DISCOVERY REACH TO NSI, 613 AND 6 

In this section we try to summarize the sensitivities to NSI, 613 and 5 that can be achieved 
by a neutrino factory with the intermediate-far detector setting. However, in looking for 



^ If one wants to interpret this symmetry as the one which exist in the truncated system in which Am^i 

is artificially switched off, it is a discrete version of the two phase degeneracy discussed in Sec. 
^ We have confirmed that the intrinsic degeneracy survives the rate only analysis. 
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ways to present the sensitivities, we recognized that the sensitivities depend very much on 
the input values of S and ^13. Therefore, we need to show the sensitivities as a function of 
both 6 and 6*13 simultaneously. 

In Figs. 18 22 we have presented the equi-uncertainty contours of a particular observable 
O in the plane spanned by the input values of sin^ 2^13 and 6. We have defined the uncer- 
tainty of measuring (or constraining) the observable O as AO = (Omax — Omin)/2, except 
for the case of sin^ 26'i3, where we give the fractional uncertainty, i.e. A(sin^ 2^13)/ sin^ 2^13. 
For all cases we present the 2a CL (2 DOF) contours. 

The sensitivity to the NSI parameters, 613 and 6 can be determined in the analysis of the 
6 combinations of NSI elements we have examined in this work: Eee — £e^J., £ee — ^er, £ee — ^tt, 
^e/i — £eT, £eti — ^tt and Eer — £tt- Therefore we need altogether 12 panels to fully present 
the sensitivities to the e's, and 12 panels present the sensitivities to the standard oscillation 
parameters 6^13 (6 panels) and 6 (6 panels). 

We show, for definiteness, the case the NSI parameters have zero (vanishingly small) 
input values. But, as we saw in Sec. VI, the results would not be so different even if we had 
taken non-zero input values for the NSI. Prior to showing the sensitivity contours we wish 
to warn the readers; Some of the structures are due to the finite grid used in our calculation, 
so the precise shape of the contours may not be reliable. 

Let us first look at Figs. 18 and 19 where we show the sensitivity to the NSI parameters. 
Roughly speaking, for the typical value sin^2^i3 ~ 10~^ the sensitivity to e.rT,£ee,£eT and 
Eef, are ~ 10-20%, 2-10%, 0.1-0.4% and 0.01-0.04%, respectively. As we have discussed in 
previous sections, the off-diagonal NSI parameters have a much more significant impact 
than that of the diagonal ones. We confirm here that in fact, they can be significantly more 
constrained by data. We observe that the sensitivity to the off-diagonal NSI parameters is 
basically not affected by the presence of another non-zero NSI contribution, whereas this 
is not true for the sensitivity to the diagonal elements. We suspect that this comportment 
may continue to be true even if more NSI parameters are switched on simultaneously. 

Some of the features of these contours can be readily understood from the bi-probability 
plots in Fig. [l] For example, from panel (b2) of Fig. 18 and panels (d2) and (fl) of Fig. 19 



we observe that the sensitivity to Eer is best at 5 = and 6 = tt and worst at 5 = 7r/2 and 
3n/2. This is exactly what one would expect from the upper right panel of Fig. [l| namely, 
the points in P — P space corresponding to 6 = n/2 (square) and 6 = 37r/2 (diamond) 
with non-zero Eer can be confused with that of the standard case without NSI effect (orange 
strip). One can also understand that this behavior does not depend on 613. 

Despite the contours for Eee have a more complicated structure, one can guess from the 
upper left panel of Fig. [ijthat the best sensitivity should be for 5 ~ — 7r/4 and worse for 
5 ~ vr. This is also confirmed by panels (al), (bl) and (cl) of Fig. 



The sensitivity to Eg^^ shows an intricate dependence on 6 and ^13 as seen in Fig. [T8](a2) 
and Figs. [I9|^dl) and (el). The sensitivity is the worst at around 5 ~ '7r/2 and ^ 3tt/2 
for 8111^2913 ~ 10"^ and at around 5 ~ tt for sin^ 2^13 ~ 10"^. At sin^ 2^13 ~ 10-^ the 
sensitivity to Eef_i has no significant dependence on 6. Let us understand the cause of such 
complicated features. First of all, the sensitivity to Eefi is coming dominantly from the 
measurement at 3000 km, as indicted in Figs. |4| |6] and 12 Then, it should be possible 



to understand the behavior at sin^ 26*13 ~ 10~^ from the bi-probability plot give in Fig|2} 
This is indeed the case; The upper middle panel of the bi-probability plot shows that the 
ellipses of positive and negative Ee^ as well as the one without NSI overlap with each other 
at 5 ~ TT, and hence the sensitivity is worst at 5 ~ vr. The similar bi-probability plot for 
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sin^ 26'i3 ~ 10~^ (not shown) indicates that a positive, negative, and zero Ee^ elhpses meet at 
5 ~ 7r/2 and ~ 371/2. (Precise vahies of the crossing point are at around 6 = O.Gvr and lAn.) 
At sin^ 2^13 ~ lO""^ (not shown) there is no particular value of 6 at which these ellipses 
overlap. Thus, these features of the bi-probability plots explain the behavior of sensitivity 
contours of ^e/i in Fig- [l8|^a2) and Figs. [T9|^dl) and (el). 

Let us also discuss the sensitivity reach of ^13 and 6 with and without NSI effect. In 
Fig. 20 we show the iso-contours of uncertainty on the determination of the CP phase 6 
(upper panel) and of sin^ 26*13 (lower panel) in the plane spanned by the input (true) values 
of sin^ 2^13 and 6 in the absence of NSI effect. These results must be compared with the 
case with NSI effect we will show below. 

In Fig. 21 we show the iso-uncertainty contour of the CP phase 6 in the presence of NSI 
parameters for the 6 different combinations of the 2 e systems we considered in this work. 
Roughly speaking, for the typical input value sin^ 2^13 ~ 10~^ the sensitive become worse 
when NSI is included. Nevertheless, the change is from A6 ^ 0.05 — 0.1 radians (without 
NSI) to ^ 0.1 — 0.15 radians or so (with NSI), in all 6 cases. We conclude that the difference 
between the sensitivities to 6 with and without NSI is not dramatic. 

On the other hand, regarding the sensitivity to sin^ 26*13, comparing the lower panel of 
Fig. 20 with Fig. 22 , we can see that for the typical input value sin^ 2^13 ~ 10"^ the fractional 
uncertainty on sin 2^13 becomes larger with NSI, going from A (sin 2^13)/ sin^ 2^13 ~ 10% 
(without NSI) to ~ 10 — 20% (with NSI). Again we conclude that the impact of NSI in the 
determination of ^13 is not so striking. 



VIII. CONCLUDING REMARKS 

We have demonstrated in this paper that a neutrino factory equipped with an intense 
neutrino flux from a muon storage ring and two detectors, one located at L = 3000 km and 
the other at L = 7000 km, is powerful enough to probe into extremely small values of the NSI 
parameters. We have relied on the golden channel, Ug — > u^, and its anti-neutrino counter 
part in our analysis in this paper. Six different combinations of the two e systems that can 
be obtained from the four NSI parameters, ea/3 with (a,/5) being (e,e), (r, r), {e,fi) and 
(e,r), are analyzed under the assumption of ignoring the effects of NSI at the production 
and the detection of neutrinos. 

The sensitivities to off-diagonal e's are excellent, l^erl — a few x 10^'^ and l^e^d — a few x 
10""^ and while the ones for the diagonal e's are acceptable, |eee|(krr|) — 0.1(0.2) at 3a CL 
and 2 DOF. These sensitivities remain more or less independent of ^13 down to extremely 
small values such as sin^ 26*13 = 10~^. They seem also very robust in the sense that they 
are not very disturbed by the presence of another non-zero NSI contribution. The above 
characteristics of the sensitivities to NSI suggest that in our setting the off-diagonal e's are 
likely the best place to discover NSI. We note that these results are obtained under the 
assumption of ignoring background as well as systematic errors in our analysis. According 
to our estimate, however, the effect of inclusion of them is limited to ~ 50% (a factor of 2) 
or so for ^13 as small as sin^ 2^13 = 10"'^ (10~^). 

One of the most significant features of the results obtained in our analysis is that the 
presence of NSI does not confuses the precision measurement of ^13 and S. The favorite 
property arises from the synergy between the two detectors placed at the two different 
baselines; The detector at the magic baseline is extremely (reasonably) sensitive to the 
off-diagonal (diagonal) e's but lack sensitivity to 6. On the other hand, the intermediate 
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detector at L = 3000 km is sensitive to S, while lacking good sensitivity to e's (except in the 
case of £en)- We have shown in Sec. |V] and VI that, when combined, the synergy between 
the two detectors has an enormous power to resolve the confusion between ^13 and NSI. 
Moreover, the impact of the systematic errors on S and ^13 is limited especially at small ^13 
and at most ~10%-20% level at sin^ 261^ = 10~^. We believe that the results obtained in 
this paper open the door to the possibility of using neutrino factory as a discovery machine 
for NSI. 

Our analysis in this paper, however, has limitations of validity because we have assumed 
that the effects of NSI in production and the detection processes are negligible. It may be 
a self-consistent approximation when we discuss the Eer system, assuming a very small Eer, 
because its effect on muon decay is expected to be small. But, if we discuss the e^^ system 
it affects the production and the detection processes at one-loop level and all these effects 
have to be consistently dealt with. In this paper we are not able to address this point. 
Our analysis also does not include the energy resolution and the uncertainty on the matter 
density as well as on the remaining mixing parameters. 

We have also mentioned in Sec.|Tl]that the phase of the off-diagonal Ea/3 produces confusion 
problems of two types, the two-phase and the phase-magnitude confusion. In our analysis we 
have shown that the discrete version of the former degeneracy is resolved by the two-detector 
setting examined in this paper. Therefore, it is natural to expect that the degeneracies in 
its generic form as well as the latter type will also be resolved by the two-detector setting. 
The magnitudes of the solar Arn^i sensitive terms will be different between the intermediate 
and the far detectors so that the confusions will be lifted by combining them. The work 
toward this direction is in progress. 

One could think about adding the silver channel in the analysis [IE]. This could, in 
principle, increase the sensitivity to e^x- In the present setting of magnetized iron detectors it 
requires separation of muons produced by tau decay from those of oscillated z/^ CC reaction. 
This requires more sophisticated analysis and therefore it is left for future investigation. 
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sin^2e,i and 5 marginalized ■"''"'= .i"X,='> »»i 



£ =E ~£ -0 
« ET ru 




J 0,1 0. 



FIG. 3: Allowed regions projected into the plane of 2 NSI parameters, Eee-^er (left panels), Err-^er 
(middle panels) and eee-^rr (right panels) corresponding to the case where the input parameters 
are sin^ 2^13 = 0.001 and 6 = ir/A and no non-standard interactions (or all the e's are zero), for 
= 50 GeV and the baseline of L = 3000 km (upper panels), 7000 km (middle horizontal panels) and 
combination (lower panels). The thin dashed lines are to indicate the input values of Eap- The fit 
was performed by varying freely 4 parameters, ^13, 6 and 2 e's with ^13 and 6 being marginalized. 
The number of muons decays per year is 10^^, the exposure considered is 4 (4) years for neutrino 
(anti-neutrino), and each detector mass is assumed to be 50 kton. The number of energy bins 

2.5x10-3 eV2, 

8.0 X 10"^ eV^ and sin^ 6*12 = 0.31. 



considered is three. The other standard oscillation parameters are fixed as Am23 



sin^ 6*23 = 0.5, Amf2 
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FIG. 4: The S3-II16 3/8 in Fig. jsjbut for a, difFcrent combination of 2 6's, s^^-s^^ (left panels), Srr-e 
(middle panels) and ee/i-Eer (right panels). 
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FIG. 6: The same as in Fig. |4]but with 6 = 37r/2. 
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FIG. 7: Allowed regions projected into the plane of sin^ 20i3-(5 corresponding to the case where 
the input parameters are sin^ 2^13 = 0.001 and 6 = tt/4 and no non-standard interactions (or 
all the e's are zero), for E"^ = 50 GeV and the baseline of L = 3000 km (upper panels), 7000 
km (middle horizontal panels) and combination (lower panels). The fit was performed by varying 
freely 4 parameters, ^13, 6 and 2 e's where See and Eer are marginalized (left panels). Err and Eer 
are marginalized (middle panels) and See and Err are marginalized (right panels). The same input 
and fitting parameters as in Fig. |3]but projections are made into the different parameter space. 
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FIG. 8: The same as in Fig. [7]but for different combination of 2 e's to which the fit to sin^ 2^13 
and 6 (left panels), ErT-^en (middle panels) and EefM-Eer (right panels). 
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FIG. 10: The same as in Fig. [s] but with 5 = 37r/2. 
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FIG. 11: The same as in Fig. [5] but with sin^ 26^13 = 0.0001. 
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FIG. 12: The same as in Fig. [g] but with sin^ 26^13 = 0.0001. 
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FIG. 14: The same as in Fig. 10 but with sin^ 29ri = 0.0001 
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NoNSI - sin^2e., = 0.001 




FIG. 15: Allowed regions in the plane sin^ 2^13 — S for the standard oscillation case (no NSI) for 
the input values sin^ 2^13 = 0.001, 6 = 7r/4 (left panels) and 6 = 37r/2 (right panels). The allowed 
regions were computed for 2 DOF. The experimental conditions are the same as for Fig. [7j 
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FIG. 16: The same as in Fig. |5]but for non-vanishing input values of e; e^e = 0.1, Sgr = 0.01 and 
= 0.2. We note that only the input values of 2 e's are set to be non-zero at the same time. 
The thin dashed lines indicate the corresponding non-zero values of for each panel. 
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FIG. 17: The same as in Fig. |9]but for non-vanishing input values of e; £ee = 0.1, Ser = 0.01 and 

Err = 0.2. 
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FIG. 18: Isocontours of 2 a CL (for 2 DOF) sensitivity (uncertainty) for 8^/3 parameter as 
functions of the input parameters of sin^ 2^13 and S. For each point in the sin^ 2^13 — S plane 
the uncertainty is defined as Ae = (£max — £^min)/2 where £inax(min) indicates the maximum and 
minimum allowed value of e parameters which is consistent with the case without NSI effect. In the 
upper, lower and bottom panels, the sensitivities for (al) Sge and (a2) for the Eee-^e^ system, 
(bl) See and (b2) for the e^e-^eT system, and (cl) e^e and (c2) Et-t for the eee-^rr system, 
respectively, are shown. We note that the the uncertainty for shown in (a2) and £er shown in 
(b2) are magnified by 100 and 10, respectively. Some of the structures are due to the finite grid 
used in our calculation, so the precise shape of the contours may not be reliable. 
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FIG. 19: Same as in Fig. 18 but for different combination of the 2 e system. In the upper, lower 
and bottom panels, the sensitivities for (dl) and (d2) Ser for the eef_i-£eT system, (el) Eefi and 
(e2) Srr for the e 



system, and (fl) Eer and (12) Err for the Eer-^rr system, respectively, are 
shown. As in Fig. 18 the uncertainty for shown in (dl) and (el) and Eer shown in (d2) and 
(fl) are magnified by 100 and 10, respectively. 



36 



5 - sensitivity 
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FIG. 20: Upper panel: Iso-contours of 2 o" CL (for 2 DOF) sensitivity (uncertainty) for the CP 
phase 6 (in radians) expected to be achieved at neutrino factory in the absence of the NSI effect in 
the plane of the true values of 6 and sin^ 2^13. The uncertainty is defined as A6 = ((5max — '^min)/2 
in radians, where (Jmax(min) is maximum (minimum) allowed value of 6 (mod. 2tt) for each given 
input point. Lower panel: Similar plot as in the upper panel but for the fractional uncertainty 
A(sin^ 2^13)/ sin^ 2^13 (in percent) is shown. 
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5 - sensitivity 




FIG. 21: Same as in the upper panel of Fig. [20] but for the case where the NSI parameters are 
turned on in the fit. The iso-contours of A6 (in radians) at 2 cj CL (for 2 DOF) are shown for the 

6 combinations of 2 e system: (a) See - £e^l, (b) See - ^er, (c) See - ^tt, (d) Se^ - Eer, (e) ee/i - Err 
and (f) Eer — Err- 
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FIG. 22: Same as in the lower panel of Fig. [20] but for the case where the NSI parameters are 
turned on in the fit. The iso- contours of A(sin2 26*13)/ sin^ 26*13 at 2 a CL (for 2 DOF) are shown 
for the 6 combinations of 2 e system: (a) See — Ee/i, (b) See — ^er, (c) £ee — £tt, (d) £e^i — ^er, (e) 

^efi ~ ^TT and (f) Eer ~ ^tt- 
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APPENDIX A: DERIVATION OF APPEARANCE PROBABILITY P{i^e ^ i^f,) 
IN THE PRESENCE OF NON-STANDARD INTERACTIONS 



In this Appendix, we give a self-contained discussion for deriving the expression of the 
appearance probabihty P(z^e ^^i) with simultaneous presence of non-standard interactions 
See and Est- For simplicity, we denote it Eer — £ee system, and use the method developed 
by Kimura, Takamura, and Yokomakura (KTY) [49J. (See [50] for minor sign error in the 
original formula, and [40] for a reformulation.) 

The evolution equation of neutrinos can be written in the flavor eigenstate as 



d 

ax 



1 

2E 



{a, 13 = e,/i,r). 



(Al) 



where the Hamiltonian is given by 



H = U 


















+ a{x) 



1 -\- 6^gg Sf^j- 



E* 



whose first term will be denoted as i/™^ hereafter and Ami- = — ^ 



(A2) 







f. (Hence, Amf^ 

by definition.) In (A2), a = 2\/2Gpne{x)E denotes the coefficient related to the index of 
refraction of neutrinos in medium [3] with electron number density n^lx), where Gp is the 
Fermi constant and E is the neutrino energy. Despite that ne{x) may depend upon locations 
along the neutrino trajectory, we use constant density approximation throughout this paper. 
The MNS matrix U relates the flavor and the vacuum mass eigenstates as 



where i runs over 1-3. We use the standard parametrization of the MNS matrix 



(A3) 



U 



C12C13 

iS 

~'5l2C23 ~ Ci2S23Si3e 
S12S23 - Ci2C23Sl3e*'' 



S12C13 Sue 

iS 

C12C23 ~ ■5i2S23Si3e S23C13 

-C12S23 - Si2C23Sl3e*'^ C23C13 



(A4) 



where Cij and Sij 



1-3) imply cos% and sin^^jj, respectively. 
By defining renormalized matter coefficient and renormalized 



a = a(l + 



(A5) 



the matter term in the Hamiltonian can be written as 

1 Eer 



El, 

Thus, the problem is reduced to the effective system with only single type of NSI, Ee 
We now define the mass eigenstate in matter z/?" by using the transformation 



(A6) 



m 
i 1 



(AT) 
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where V is the unitary matrix which diagonahze the Hamiltonian with scaled eigenvalues A 
as V^HV = -ffdiag = diag{a\i,a\2,a\2,). We first obtain the expressions of the eigenvalues 
of the Hamiltonian (A2). They are determined by the equation det[H — Xal] = which is 
the cubic equation for the scaled eigenvalue A: 



+ 



A^ - (1 + ^31 + 52l)A' 
Clsl^Sl + {^31 + Ci2 + + 2Ci2Si2S23Ci3i?e(£er)} (^21 

-2C23Ci3Si3(53i - sl2S2l)Re{eere'^) - l^erP 



+ e 



2 



C12C13 + 2Ci2Si2S23Cl3Re{6er) - 2Ci2C23Cl3Sl3-Re(ee^e*' 



0. 



(A8) 



'S23Cl3'^31 + (C12C23 + ■5l2'523'5l3 " 2C12S12C23S23S13 COS (5)521 

in which everything is scaled by a and 621 and 531 denote the scaled squared mass differences, 

021 = — z — , ^31 = — z — • (A9) 



1. KTY method for obtaining exact oscillation probability with NSI 

We follow the KTY method [IH] for deriving P(t'e ~^ '^n) and write down the equations 



TT TTvac 

TT TT U U ( Zjuac I ~ N TTvac TTvac Tjvac 

They give relationships between mixing matrix in vacuum and in matter as 

i i 
cyclic cyclic 

J2 >^j^kVeiV;, = J2 Sji^kiUeiU;, + J2 ^^^UriU;, = q. 

ijk ijk i 



(AlO) 



(All) 



Notice that the effect of Eer is contained only in q. Solving (All) for KjV^* under the 
constraint of unitarity ^ • VeiV*^ = we obtain 



v.v*. = Ph±l 



(A12) 



where Aji = \j — Aj and k) are cyclic. 

Then, the appearance probability P(z/e Ue) is given exactly by 



cyclic 



aL 



4 Yl (^Jel + Re Jii) cos [ — A,, j sin [ — A,, j sin — A,, 

(ijk) 
cyclic 



aL 



aL 



(ijk) 



aL 



E J^'"" ( j sin f ^A23 ) sin ( ^A3i ) , (A13) 



. / aL 
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where the sum over the cychc permutation is imphed and 



J 



lm{pq* 



A12A23A31 



(A14) 
(A15) 



It is thus convenient to compute the combination Re J*^ + Re J^^ (note the minus sign 



e/i 



J. 



(A16) 



where 



Jj = \p\'X'j + 2Reipq*)Xj + \q\ 



Ci + 2Ai^^ cos 5 + 2 AY'' cos 26 + 2B\'' sin 5 + 2B\''> sin 25 



in. 



Cj = (Po+K)\ +2{poRe(go) +PiRe(gi)}Aj + |go| + |gi| + |g2| 



(A17) 



A 
A 



B 



(I) 
j 

in) 
j 

in 



ill) 



PoPi^j + {PoRe(gi + q2) + PiRe(go)}Aj + Re(go)Re(gi + ^2) + Im(go)Im(gi + ga) 
PiRe(g2)Aj + Re(gi)Re(g2) + Im(gi)Im(g2) 

{poIm(gi - ^2) - piIm(go)}Aj + Re(go)Im(gi - ^2) - Im(go)Re(gi - ^2) 
-piIm(g2)Aj + Im(gi)Re(g2) - Re(gi)Im(g2) 



J is given by 
J 

where 



1 



A12A23A31 



J(^) sin 6 + J(^^) sin 25 + iT^^^ + K^^'^ cos 5 + iT^^^^ cos 25 



(A18) 



(A19) 



i^(o) 



{poRe(gi - 52) - PiRe(go)} 
-piRe(g2) 
-polmgo - Pilmgi 
-{poIni(gi + ^2) +PiIm(go)} 
-pilm(g2) 



The coefficients p and g, which are defined in (All), can be written as 

i5 



P = Po+ Pie 



go + qie-'' + q2e+'' 



(A20) 



(A21) 



where 



PO = 521C12S12C23C13, 

Pi = {S3I - Si252l)s23Cl3Sl3 

-<531'^2lCl2Sl2C23Cl3 + ^erC23S23 



Ql 
12 



^3lCi3 - 52l(Ci2 - S'oS 



2 2 N 
12*13/' 



5. 



21 



'^3lCi2S23Cl3Sl3 + £^eTCl2Si2S23Si3 



~ r 2 
-£eTf2lCl2'Si2C23Si3. 



(A22) 
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Notice that the coefficient p is identical with the standard case without NSI interactions 
whereas q has a httle more complex 5 dependence with q2 term and the coefficients qi 
= 1 — 3) have imaginary parts. 



Collecting formulae given in the equations from (A13) to (A22) and using the exact eigen- 



values by solving the cubic equation (A8), we obtain the exact expression of the appearance 
probability P{ve —>■ v^) with the neutrino NSI Eee and Eer- Notice that a and Eer are the 



renormalized quantities defined by (A5) 



cyclic ^ 

= 8 y ~ \-C, + Af cos 5 + cos 25 + ^ sin 5 + sf^ 
' aL . \ . f dL . \ . f dL 



sin 25 



, aL \ . f aL \ . f aL \ 
^ (^Ie'^'^ ) ^'"^ [lE^'^J '''' (,4E^^^ ) ' 

. ^ [jW sin5 + J(^^) sin 25 + K^'^ + K^'^ 

^12^23^31 - 

X sin ( — Aio I sin ( - — A23 | sin | ^— Aqi | . 

\AE \AE V \4:E ^ 



+ 8 



cos cos 2(5 



(A23) 



2. Leading order formula of P^i^e '^n) 

Now, we define the perturbative scheme we use in our derivation of the approximate 
expression of the appearance probability P(z/e J-'^)- We regard ^si and a as of order 
unity, and assume that S13 ~ 821 — Eer are small and are of the same order to organize the 
perturbation theory. We denote the expansion parameters symbolically as e. We assume 
that e ~ 10~^. If we organize the perturbation expansion in terms of e we will recognize 
that the leading order terms in the oscillation probability in P{Vf. u^) is of order e^. 

To the next to leading order in e the solutions of the equation are given under the 
convention that Ai < A2 < A3 by 

A2 = 531--^|4 + 2C23Sl3Re(£ere*')| 
1 - ()31 L J 

A3 = 1 + Y^,{^^^^ + 2c23Si3Re(£ere^')} + SI2S21 (A24) 



In (A24) we have ignored even smaller corrections to the lowest eigenvalue because it is of 
order"?. 

We restrict ourselves to the leading order as was done by Cervera et al. J28j for the 
standard case without ^o/j's. To this order, we in fact do not need the order e corrections 
in (A24). By putting back the physical quantities replacing the scaled variables, the z/^ 



43 



appearance probability can be expressed to order as 
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(A25) 



where P{i'e ^^1]^ = ^)\2nd is nothing but the Cervera et al. formula 
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cos (5-^ 



(A26) 



It is notable that the NSI effects survive in the leading order, ~ e^. If one want to have 
explicit expression with e^t and See one can just use the relations (A5) in (A25). The formula 
for P{ve i^fi) is valid only for small Est but for any finite size See- The antineutrino 
probability can be obtained by the replacement 6 — > —6, a —a, and Ea/s — >■ e^js- We have 
checked that the same analytic formulas are obtained, when expressed in terms of observable 
physical quantities, even if we work in the intermediate energy region where A2 > A3. 
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Similarly, the formula of P(z/e t'^) with Sen can be computed to the leading order as 
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